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Abstract 

We describe a method for generating entanglement between two spatially separated dipoles 
coupled to optical micro-cavities. The protocol works even when the dipoles have different resonant 
frequencies and radiative lifetimes. This method is particularly important for solid-state emitters, 
such as quantum dots, which suffer from large inhomogeneous broadening. We show that high 
fidelities can be obtained over a large dipole detuning range without significant loss of efficiency. 
We analyze the impact of higher order photon number states and cavity resonance mismatch on 
the performance of the protocol. 

PACS numbers: 



I. 1. INTRODUCTION 



Generation of entanglement between qubits is an important operation for a large variety 
of applications in quantum information processing. Such states can be used to the real- 
ize schemes such as transmission of secret messages via quantum key distribution P, [2] and 
teleportation of quantum information [3j HI El E]. The exchange of entanglement between 
two distant parties is also required for implementation of quantum repeaters |7J which use a 
combination of entanglement swapping and entanglement purification [8 J to achieve uncon- 
ditional secure communication over arbitrarily long distances. 

To date, a variety of methods have been proposed for creating entanglement between 
spatially separated nodes. One of the most common methods is to transmit entangled 
photons generated by parametric down-conversion [9]. Entanglement protocols for atomic 
systems have also been proposed [10J HU [121 021 EH]- Atom entanglement has the advantages 
that quantum information can be stored for long time periods, which is important for long 
distance quantum networking. 

Semiconductor based approaches to quantum information processing are currently an 
area of great interest because they offer the potential for a compact and scalable quantum 
information architecture. Furthermore, solid-state emitters such as semiconductor quantum 
dots (QDs), can be coupled to ultra-compact cavity waveguide systems to form highly inte- 
grated quantum systems US] - A major challenge in using solid-state emitters is that they 
suffer from enormous inhomogeneous broadening, typically caused by emitter size variation 
and strain fields in the host material. The inhomogeneous broadening makes it difficult to 
find two emitters with identical emission wavelengths. Protocols to date for generating atom 
entanglement require the dipoles to emit indistinguishable photons, and are thus difficult 
to implement in semiconductor systems. In order to implement quantum networking in 
semiconductors, we need a protocol that works even when the dipoles emit photons that are 
distinguishable. 

In this paper, we describe a protocol for creating entanglement between two dipoles with 
different radiative properties, such as different emission wavelengths or radiative lifetimes. 
The proposed protocol uses Dipole Induced Transparency (DIT) to achieve the desired entan- 
glement which occurs when a dipole is coupled to an optical cavity [17]. When the coupling 
is sufficiently strong, the dipole can switch a cavity from being highly transmitting to highly 
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reflecting. The switching contrast is determined by the atomic cooperativity, which is the 
ratio of the lifetime of the uncoupled emitter to the modified lifetime of the cavity-coupled 
emitter. Enhancement of spontaneous emission has been observed in semiconductor emitters 
coupled to a variety of different micro-cavity architectures [T6l [18], [T9l |20| EQ G2]. Modifica- 
tion of cavity reflectivity by coupling a quantum dot to a photonic crystal nanocavity has 
been recently observed \23\ [21]. 

In section 2, we describe the basic protocol under idealized assumptions that all fields can 
be expanded to first order in photon number and that the two cavities have the same reso- 
nant frequencies. Section 3 then considers the effect of higher order photon number states 
on the efficiency and fidelity of entanglement. The impact of non-linear behavior away from 
the weak excitation limit is also discussed. In section 4 we investigate the effect of cavity 
frequency mismatch on the entanglement. Finally, in section 5 we perform a precise numeri- 
cal simulation of the entanglement generation protocol for one specific implementation using 
the exciton bi-exciton cascade of a single Indium Arsenide (InAs) quantum dot. Numerical 
results from recent experimental work is used to show that this may be a promising method 
for achieving entanglement between two QDs for the first time. 

II. 2. PROTOCOL FOR ENTANGLEMENT GENERATION 

QD Entanglement 




FIG. 1: Schematic of cavity waveguide system for generating entanglement between two spatially 
separated dipoles using DIT 

The schematic for generating entanglement between two spatially separated dipoles that 
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emit distinguishable photons is shown in Fig 1. Each qubit consists of a dipole coupled to 
a double sided cavity. Each dipole is assumed to have three states: a ground state, a long 
lived metastable state and an excited state, which we refer to as \g), \m) and |e) respectively. 
The states \g) and \m) represent the two qubit states of the dipole. The transition from 
the ground state to the excited state for dipole 1, may be detuned by 8\ from the resonant 
frequency uj\ of cavity 1. Similarly, the transition from the ground state to the excited state 
for dipole 2 may be detuned by 82 from the resonant frequency 002 of cavity 2. We assume 
that when the dipole is in state \g) it couples to the cavity mode via an optical transition 
to state |e), while when it is in state \m) it does not optically couple to the cavity mode. 
State \m) may be decoupled from the cavity due to either spectral detuning or selection 
rules. Although state \m) is illustrated in the diagram as having an energy level that is 
in between states \g) and \e), this is not required. The only requirement is that when the 
dipole is in state \m) it is decoupled from the cavity. This point will be analyzed in more 
detail when we consider using the exciton-biexciton transitions of an InAs QD to achieve 
entanglement. The desired level structure described above can be realized in a variety 
of solid-state material systems. In semiconductor quantum dots one can use the exciton 
and biexciton transitions [25], as well as the spin-based bright and dark exciton states[26j. 
In addition, three level structures can also be achieved using quantum-dot molecules pT], 
charged quantum dots[28] and impurity bound excitons[29]. Similar qubits states could also 
be realized in other materials such as diamond using neutral and negatively charged nitrogen 
vacancy defects [SHI EI]- 

The decay rates of the two dipoles is given by 71 and 72 respectively. To characterize the 
interaction between the dipoles and the cavity modes, we define the operators <t£_ and 02- ■ 
They represent the dipole lowering operators for the dipoles in cavities 1 and 2 respectively. 
It should be noted that o"i_ and 02- represent the dipole lowering operators for the g- 
e transition. Although state \m) is decoupled from the cavity, we still define the dipole 
operators cr m i_ and (7 m V for the dipole in state \m). These dipole are detuned by 8 m \ and 
8 m 2 from their respective cavities. 

We define a in and c in as the two input modes, a oni and c out as the reflected modes, and 
h out and d out as the transmitted modes to the two cavities, as illustrated in Fig. 1. The 
energy decay rate of cavity 1 into the reflected and transmitted modes is given by K r \ and K t i 
respectively. Similarly, the energy decay rates of cavity 2 into the reflected and transmitted 
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modes is given by K r 2 and n t 2 respectively. There is also the decay rate kh and K12 into 
the parasitic leaky modes that is due to losses such as material absorption and out of plane 
scattering. The field inside the cavities are represented by the cavity field operators fi and 
f 2 - 

The protocol works as follows. Both the dipoles are initialized to be in an equal super- 
position of qubit states \g) and \m). This can be achieved by first driving the dipoles into 
the lowest energy state by either waiting several radiative lifetimes or optical pumping. The 
qubit state can then be rotated by either a direct 7r/2 transition, or a raman transition[32j. 
The choice of method depends on the specifics of the dipole and material system. Once 
the initialization step is complete, the initialized state of the two dipole system is given by 
1/2(|00) + \mm) + \gm) + \mg)). 

After the initialization of the dipoles, a weak coherent field \a) with frequency u is inserted 
at input k in . Simultaneously, another weak coherent field \/3) that is phase coherent with \a) 
(i.e. originates from a common laser source) is injected at c in . These input fields interact 
with the cavity-dipole system. The interaction between the input field a iri and cavity-dipole 
system 1 can be characterized by the Heisenbergs equations of motion for the cavity field 
operator fi and the dipole lowering operator <j\- 



— = -(zo) + ( K n + «ti + «n)/2)fi - v^i^m - ig&i- ~ W^tal- 
{ ' ( ' -- (-*( w o + <*i) + 7i)o"i- + iga A fi (!) 



dt 



-i(u + 5 m l) + 7i)cr m i_ + iga zml ft 



dt 

Similar equations can also be written for the interaction of the input field Cj n with cavity- 
dipole system 2. 

The interaction between the input fields and the cavity-dipole systems results in part of 
the field being transmitted into the modes Y) out and d^ ui , while the remainder is reflected 
into the modes a^ Mi and cl ut , or absorbed by the QD. The amount of light reflected and 
transmitted is given by the cavity reflection and transmission coefficients. Our analysis 
works in the weak excitation limit, where predominantly the quantum dots are populated 
in the ground state. In this limit, (<r z i(t)) ~ —1. This also implies that the population 
inversion for the dipole in state \m) is close to i.e. (cr zm i(t)) ~ Using this limit in Eq. 1, 
we can derive the reflection and transmission coefficients to be [331 



5 



, v {-^1 + ^(A^ ) + («rl ~ ^1 ~ «ll)/2 

= — 2 

(-iAwx + (« rt + Kti + k,i)/2 + _ <(2Wl £ fl|3)+71 ) (2) 

tl(Cj) = 2 

(-iAu! + (/en + K tl + «n)/2 + _ t(AcJl 3 „5 l)+7l ) 

where Ac^i = lo—uj\. These equations are obtained for cavity-dipole system 1. The reflection 
and transmission coefficients for dipole 2 are identical to dipole 1 in form, and are obtained 
by substituting the dipole 2 parameters into Eq. 2. 

To get a better feel for Eq. 2, it is helpful to first consider the simplified case where 
Aui = 0. Assume first that the dipole is in state \g) and that the g-e transition is resonant 
with the cavity such that Si = 0. We see that maximum reflection and minimum transmission 
occurs for the case when K r i = K t i+Kn, called the critical coupling condition. This condition 
ensures that no light is reflected from the cavity when the incident field is directly on cavity 
resonance. We represent the decay rate K r i at critical coupling as K\. Hence, the transmission 
and reflection coefficients simplify to t = 1/(1+C) and r = C/(l+C), where C=g 2 /71/C1 is 
called the atomic cooperativity. If C>> 1, which is the desired operation regime, then r=l 
and all of the light is reflected. Now suppose the dipole is instead in state \m) which is 
detuned from the cavity by S m . We then have t = 1/(1+CL) and r = CL/(1+CL), where 
L=7i/(7i + iS m ) is a Lorentzian function. If we assume that either state \m) is highly 
detuned from the resonance of the g-e transition (S m » g 2 7 ' Ki) or the transition is very 
weak due to selection rules (C~ 0), then t=l and now the light is completely transmitted. 
Thus, by changing the state of the dipole from \g) to \m) we can completely change the 
reflectivity of the cavity. 

In a realistic system we cannot assume that the two dipoles are resonant with their 
respective cavities, since in general they will have different resonant frequencies. Nor can we 
assume the reflection and transmission coefficients will reach their ideal limits because S m is 
not infinitely large and we usually don't have perfect selection rules to cancel out the m-e 
transition. In this case we define for dipole 1, tf and r^ as the transmission and reflection 
coefficients when the dipole is in state \g), and t™ and r™ when the dipole is in state \m). 
We define rf ,tf ,r™, and t™ analogously for dipole 2. These coefficients can be calculated by 
plugging in the appropriate values corresponding to the different transitions of the dipoles. 
It is important to emphasize that we do not assume that the S, g, 7, and k are the same for 
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both dipoles. The protocol we describe works even if all of these parameters are different, 
which is why it is so useful in semiconductors. 

Before continuing, it is worth noting that in much of the literature the atomic coopera- 
tivity C is often interchanged with the Purcell factor, defined as the ratio of the lifetime of 
the dipole inside the cavity to that of a dipole in bulk or free space, which we denote 'jbuik- 
Although the atomic cooperativity is the correct parameter to use in the strictest sense, it 
is a very difficult parameter to measure. The Purcell factor in contrast is easier to measure 
and almost always gives a lower bound on the parameter C in any realistic system. The 
reason for this is that % u ik is due both to radiative and non-radiative decay. In contrast 
the decay rate 7 is the decay rate into non-cavity modes and is mainly due to non-radiative 
processes, as radiation into modes other than the cavity mode is highly suppressed. Thus, 
outside of some atypical cases where the cavity has two modes resonant with the QD, we 
expect that ^buik > 7- Therefore, in virtually all cases C can be replaced by the Purcell 
factor to get a lower bound on the performance of the system. 

We first investigate the protocol under the assumptions that the resonant frequencies of 
both the cavities are the same(wi = U2), and the input fields \a) and \/3) are sufficiently 
weak that we may expand them to first order in photon number. The initial state of the 
system(dipoles and fields) is given by j^) = l/2(\gg) + \mm) + \gm) + \mg))(aa. in + f3c in ). 
The fields, after interacting with the cavities, are transformed according to cavity reflection 
and transmission coefficients. That is, if dipole 1 is in state \g) then k in — > rfk out + t{b out , 
and if it is in state \m) then a in — * r™k out + t™b oui . The transformation for photon in 
c in and dipole 2 is defined in a completely analogous way. The reflected fields from the 
two cavities are mixed on a 50/50 beamsplitter that applies the transformation: — > 
(di + d 2 ) /V%, c\ ut -> (di - d 2 ) /V2- The final state of the QDs can be obtained by applying 
the cavity and beamsplitter transformations. If a detection event is observed in detector d 2 , 
then the state of the two QDs collapses to 



l*/> = ^(K? - fri)\99) + ur{\gm) - [3r 9 2 \mg)] (3) 
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where N 2 = \ar{ — I3r g 2 \ 2 + \arf\ 2 + \(3r\ 

In general, rf 7^ rf because the dipoles have different resonant frequencies. However, 
we can correct for this mismatch by properly adjusting the amplitudes of the fields. If the 



amplitudes are selected such that 

arf = (34 (4) 

the state of the qubits is projected onto |\&_) = (\gm) — \mg))J v2 which is an ideal entangled 
state. Thus, by properly choosing the amplitude and phase of the input coherent fields \a) 
and \/3), we ensure that a detection at d 2 creates an entangled state of dipoles. Note that 
the entanglement generation is accomplished despite the fact that the two dipoles may have 
completely different resonant frequencies or decay rates. 

III. 3. HIGHER ORDER PHOTON NUMBERS 

The matching condition ar\ = (5r\ as given in Eq [4] gives a relationship between a and 
(3, but does not tell us how large to make a. In general, we want to make \a\ 2 as large as 
possible to improve the chances of a detection event at d 2 . The probability of detecting 
a photon at detector d 2 is defined as the efficiency rj of the protocol. When the fields \a) 
and \(3) are weak, efficiency of the protocol is proportional to the intensity of the field at 
d~2 and can be derived to be |ar^| 2 /4. The factor 1/4 appears because 50% of the field is 
transmitted into the modes h? out and d^ Mt and another 50% is lost when the beamsplitter 
splits the photons equally between di and d 2 . We see that we can achieve higher efficiencies 
by increasing input photon flux rate |a| 2 . However, if we make a too large we can no 
longer expand the fields to first order in photon number and higher order photon number 
contributions will become important. 

Higher order photon number contributions are undesirable because they serve as a de- 
coherence mechanism. In the ideal case where only one photon is injected into the system, 
a detection event at d 2 ensures that there are no other photons in the system which may 
carry "which path" information about the state of the dipole. Now suppose we consider the 
second order process of simultaneously injecting two photons into the input ports \a) and 
\(3). In the ideal case (both dipoles are on resonance with the cavities), if the state of the 
two dipoles is \gm), cavity 1 will reflect its incident photon while cavity 2 will transmit the 
second photon. The transmitted photon in cavity 2 will always keep track of the fact that 
dipole 2 was in state |m), and this information cannot be erased by the beamsplitter. Thus, 
we expect the state to be completely decohered when this happens. 

We will now consider not only this specific case, but full expansion of the coherent fields 
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a and (3 to all photon numbers to see how the final state of the system is affected. The 
initial state of the system is given by = l/2(\gg) + \mm) + \gm) + \mg)\a)\f3) . The 
coherent states \a) and \(3) can also be written as \ct) = Di(ct)\0) and \f3) = D 2 ((3)\0). D x 
and D2 are the displacement operators and are given by 

D 1 (a) = e *^- 01 **" 
D 2 (P) = <M,-/^» 

The displacement operator provides as convenient way of writing the coherent states and 
includes all the higher order photon numbers contributions. 

The final state of system can be obtained by applying the cavity and beamsplitter 
transformations to the initial state After applying the transformations, the final state 
of the QDs is obtained by tracing out over the photon fields conditioned on a detection event 
at detector d 2 . The state of the dipoles is therefore given by the reduced density matrix 

tr (mds) {(M\* f }(* f \M}} 
PdWOleS tr (dipolcs&ficlds) {{M\ty f )(tyf\M)} W 

The matrix M = Yl^=i \ n )d 2 ( n \ is a positive projector that projects the state of the system 
onto a subspace containing at least one photon in d 2 . This projection models the measure- 
ment performed by the photon counter, which registers a detection event as long as there is 
at least one photon in the detection mode. 

Since the final state of QDs is mixed, we need a figure of merit to measure how well the 
QDs are entangled. In this paper we use the fidelity, which is defined as the overlap integral 
between the desired final state and the actual final state of the system. In our protocol, 
the desired final state is the maximally entangled Bell state |^_). Thus, the expression for 
fidelity is (\& -\pdi P oies\^ -) • If the actual final state is same as the desired final state, we 
have a perfect entangled state and a fidelity of I. A fidelity of 0.5 implies that the state of 
the QDs is a random mixture of \gm) and \mg) and completely decohered. An analytical 
expression for the fidelity can be calculated by evaluating pdipoles and averaging over the 
state \^-)- We have carried out this calculation, but the expression for the fidelity is messy 
and the math is involved. The procedure for calculating the fidelity along with the final 
analytical expression are given in Appendix A. The expression in the appendix is used for 
subsequent calculations of fidelity. 

We also define the efficiency rj as the probability of getting a detection at detector d 2 . 
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Mathematically this is given by the expression, r\ = tr( dipo i es & fie i ds ){(M|\l//)(\l//|M)}. Using 
the matching condition arf = /3rf, the expression can be simplified to rj — 0.5(1 — e' ar i' 2 / 2 ). 

For the calculations in this paper, we use parameters that are appropriate for InAs 
quantum dots coupled to photonic crystal defect cavities. We represent the total decay rate 
out of each cavity K r \ + Kn + Kn and k t i + n t 2 + K12 as k and set it to be equal to 100 
GHz. This corresponds to a cavity Q of 3300. We set g = 20 GHz for both the quantum 
dots. We estimate dipole decay rate 7 within the cavity to be 0.125 GHz, using reported 
data that measured the lifetime of several quantum dots that were placed inside a photonic 
crystal cavity, but heavily detuned [31] . Using these values we calculate C to be 32 and the 
cavity-dipole systems to be 96.7% reflective on resonance. For the chosen values of g and k, 
the cavity-dipole systems are in the weak coupling regime(g< k/A). However, the analysis 
in this paper is completely general and is equally valid also for the strong coupling regime. 
In Fig[2j we plot both fidelity and efficiency as a function of |arf| 2 . Fidelity is plotted for 
four values of Si/k, ranging from to 1, with <5 2 fixed at 0. Note that the efficiency is only a 
function of \arf\, so the plot of efficiency is the same for all values of 81. From Fig[2j we see 
that there is a tradeoff between fidelity and efficiency as we increase a. When \arf | 2 << 1 
the fidelity is close to 1, indicating an ideal entangled state, which is consistent with our 
predictions in the weak field limit. In the region 0.1< \arf | 2 <1, the fidelity quickly drops 
due to the presence of higher photon number contributions. In the limit \arf | 2 >> 1, the 
fidelity asymptotically approaches 0.5, indicating the higher photon number contributions 
have completely decohered the state. 

When |ar^| 2 << 1, the fidelity curves for different values of 81 nearly overlap. Fidelity 
stays close to 1 in this region. However, in the region 0.1< \a\ 2 <1, the fidelity curves for 
different values of 8\ separate out. There is a drop in fidelity with increase in dipole detuning 
from Si = to 8\ = k. Also, efficiency is a function of \ar\\ alone and does not change 
with 8\. This implies that fidelity decreases with increase in dipole detuning for a constant 
efficiency. 

In Fig[2j we also plot a line of constant fidelity of 0.85. Note that for every value of 8\ 
there is a unique point on the plot corresponding to a fidelity of 0.85. As 8\ increases, this 
point shifts to lower values of |arf| 2 . Since, efficiency is a function of |arf|, this in turn 
implies a decrease in efficiency. Thus, it is important to consider how the efficiency of the 
protocol changes for a fixed value of fidelity. 
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FIG. 2: Variation of fidelity and efficiency with \arf | 2 for different values of S\. 

To investigate this, we plot efficiency as a function of 5\/k for several values of 82 for 
a constant fidelity of 0.85 in Fig [3] We see that even though there is a loss of efficiency, 
the change is gradual and there is only a 50% reduction over a cavity linewidth. Also, we 
would expect that if we added another detuning 62, efficiency would decrease. However, this 
does not happen. From Fig [3] we see that the effect of 62 is to shift the efficiency curves 
by the detuning 5 2 without altering the shape. So, the protocol can be used to obtain high 
efficiencies over a wide range of dipole detunings. 
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FIG. 3: Efficiency as a function of Si/k for different 62- Fidelity is fixed at 0.85 



85% (idelity 





IV. 4. VALIDITY OF WEAK EXCITATION LIMIT 

In the protocol we describe, dipole detuning is compensated by adjusting the amplitude 
and phase of the input coherent fields until the matching condition arf = f3r^ is satisfied. 
The more detuning we have, the larger the amplitude required by the coherent field in order 
to achieve the desired efficiency. It is possible that at some point, the amplitudes required 
by the coherent fields will be so large that the g-e transition of the QDs will be saturated 
leading to an optical nonlinearity and linewidth broadening (35]. Because of this, the cavity 
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reflection and transmission equations will depend on the pump power and Eq 2 needs to 
be modified accordingly. However, our protocol is intended to work in the linear regime 
wherein the QDs are unsaturated. This is possible only if the amplitude of the input fields 
is within a certain limit called the weak excitation limit. The weak excitation limit is defined 
as (<j z (t)) w —1, which is equivalent to the statement (er + <7_) << 1. and is necessary for 
Eq. 2 to be valid. This condition puts a constraint on the operation of the protocol. 

In order to investigate the implication of the weak excitation constraint, we start with 
the Heisenbergs equation of motion for the cavity field operator i\ and the dipole lowering 
operator a-^— given in Eq 1. We will consider cavity-dipole system 1. Similar equations are 
also applicable for cavity-dipole system 2. 

Eliminating b from Eq 1, we have 

+ 7) - 9 2 ¥- = -igV^L ( 7 ) 

Using the fact the cooperativity index C is g 2 /jk > 1, the equation can be further simplified 
and multiplied with its conjugate to obtain 

= ( ^ + \v /4) (C^) (8) 

The parameter (a~ + o~_) represents the probability of the QD being in the excited state. 
In the weak excitation limit, (cr + o~_) << 1. We also identify (aj n aj n ) as the total flux of 
photons in the input field \a). Using this in Eq|8j the weak excitation constraint thus puts 
a limit on \a\ 2 given by 

let I 2 q 2 kS 2 

— << — I h (9) 

r p k g l 

where t p is the pulse width of the laser. 

From Eq|9j we see that when there is no detuning Si, the flux of photons in the input 

2 

field | a) should be less than the modified lifetime of the QD within the cavity — . This is 
understandable because, if the first photon excites the QD and the second photon comes 
in before the QD has decayed, we will no longer be in the weak excitation limit. However, 
if the QD is off resonant from the cavity with detuning Si, not all the light that comes in 
couples to the QD. Therefore, we will be able to pump the QDs with much more power 
before we exceed the weak excitation limit. This is given by the detuning dependent term 
^in Eqg 
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Eq [9] conveys more than the weak excitation limit of \a). If we apply the matching 
condition ar\ = fir 9 2 in Eq[9j we obtain a limit on the flux of photons in the input field |/3) 
given by 

do) 

r p k g 2 

We recognize this as the weak excitation limit equation for the field \ f3) which we would have 
obtained had we used the Heisenbergs equations of motion for cavity-dipole system 2. This 
implies that if we pick \a) such that it satisfies the weak excitation limit of cavity-dipole 
system 1, the matching condition automatically ensures that the flux of photons in |/3) is 
within the weak excitation limit of cavity-dipole 2. 

Note that by making r p sufficiently long, we can always ensure that the system is in the 
weak excitation limit and that nonlinearities do no contribute. However, because we are 
using longer pulses the entanglement rate is reduced. The rate of entanglement generation 
is proportional to the rate at which the cavity reflects photons, given by R = |arf| 2 /r p . 
Using the upper limit on |a| 2 /r p from Eq[9]and cavity reflectivity rf from Eq 2, we get 

R« y - (11) 

The above equation implies that the system will remain in the linear weak excitation limit 
provided that the rate of reflected photons is less the 1 photon per modified lifetime of the 
dipole. Note that this result is true regardless of the detunings, and is therefore valid in all 
cases. 

It is instructive to compare the limits on the entanglement rate imposed by nonlinearities 
to the limits imposed by which-path information given in Section 3. The analysis of higher 
order photon numbers in the previous section showed that reflected photons \arf\ << 1 
to have a high fidelity entangled state between the QDs. In contrast, the analysis of weak 
excitation limit in this section puts an upper bound on the rate of the input photons in \a) 
and \f3) given by |arf| 2 /r p << g 2 / 'k . Thus, the two analyses are fundamentally different 
in that one limits the total number of input photons and the other limits the rate of in- 
coming photons. Although one might expect the nonlinear limit analyzed in section 4 to 
be important, it turns out that the analysis of section 3 is more restrictive, and is therefore 
the important limit to consider. To understand why, we first note that nonlinearities can 
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always be suppressed by increasing the pulse duration t p . No mater how many photons we 
inject into the system, if we make the pulses sufficiently long we will always be in the weak 
excitation limit. In contrast, which-path information does not depend on pulse duration, 
and therefore cannot be suppressed. Furthermore, in order to stay in the monochromatic 
limit (i.e. to use the single frequency approximation) it has been shown in previous work 
that the pulse duration must be longer than the modified spontaneous emission lifetime of 
the dipole[T7]. If we combine this with the results of section 3, which state that the number 
of reflected photons \arf | 2 << 1, these two conditions already constrain us to work in the 
regime where the \arf\ 2 /r p << g 2 /n. Thus, we expect the entanglement to decohere due to 
which-path information before the nonlinear behavior in section 4 is observed. For this rea- 
son deviation from weak excitation does not pose any additional restrictions to the protocol 
that were not already present in the linear scattering regime. 

V. 5. EFFECTS OF CAVITY DETUNING 

In previous sections, we considered the idealized case where both the cavities had identical 
resonant frequencies. However in realistic systems, this will not be the case. Fabrication 
imperfections may lead to slightly different resonances for the two cavities. Clearly, if even 
a small amount of mismatch between the cavities were to result in no entanglement, the 
usefulness of our protocol would be questionable. Thus, it is important to consider how 
sensitive the protocol is to cavity resonance mismatch. 

Now let's consider the case where the two cavities do not have the same resonant fre- 
quency. The analysis of the protocol in the presence of cavity detuning becomes involved 
for two reasons. First, it is no longer clear which frequency we should use for the coherent 
fields | a) and \/3). We do not know whether to place it on resonance with one of the cavities 
or somewhere in between. This can depend on both the cavity separation Au s and dipole 
detunings Si and 82. 

Second, the matching condition used in the previous section arf = fir?,, is not guaranteed 
to be optimal. If a detection event is observed in detector d 2 , then the state of the two QDs 
is 

|*/>d*ri« = ^[(orf - fri)\gg) + (ar? - f3r?)\mm) 

N (12) 

+ ar™\mg) - j3r^\gm)\ 
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where N 2 = \oer™ — f3r™\ 2 + |ar™| 2 + \ j3r™\ 2 . The matching condition arf = /3rf ensures 
that we do not have any detection at d 2 if both the dipoles are in the state \g). However, 
the field amplitude at d 2 if both the dipoles are in the state \m) i.e. (ar™ — fir™) is not 
compensated. This results in imperfect destructive interference at detector d 2 . Thus, there 
is a small probability of detection at d 2 when both the dipoles are in state \m). This causes 
a loss of fidelity. In order to obtain the state that comes closest to the desired final entangled 
state, we must optimize the fidelity with respect to u, a and (3. 

For calculating the effects of cavity detuning, we choose the frequency midway between 
the two cavity frequencies as the reference frequency A re f . Based on this reference frequency, 
uj\ = —Auj s /2 and co> 2 = Auj s /2. Also, it will be easier if we define the dipole detunings in 
terms of the reference frequency rather than the cavity frequencies. We define Ai = 5± + uo\ 
and A 2 = 62 + 002, which are the dipole detunings of dipoles 1 and 2 with respect to the 
reference frequency located midway between the two cavities. These definitions ensure that 
when increasing the cavity separation Auj s we do not affect the QDs. This is important 
because we can obtain information about the effects of cavity detuning alone by making 
these definitions. 

Figure 4 plots the dependence of fidelity on the laser frequency for several different values 
of Ai. The cavity separation Au s = 002 — oj\ is set to 50 GHz, and A 2 = 0.25k. The figure is 
optimized over the real and imaginary parts of %. The value of the maximum fidelity for the 
three curves occurs at three different frequencies. The frequency at which we get maximum 
fidelity is the optimal frequency uj. The fidelity at that frequency is the maximum fidelity 
that can be obtained for that particular configuration of Au s , A\ and A 2 . 




Aco/k 

FIG. 4: a) Fidelity as a function of laser frequency for different values of Ai. Optimization is 
performed over the real the imaginary parts of §. Cavity separation Aw s is set to 50 GHz and 
A 2 = 0.25k GHz 



15 



In Fig. 5, we plot optimized fidelity as a function of cavity detuning Au) s for different 
values of Ai with A 2 = 0. When Auj s = 0, which represents the case when there is no cavity 
detuning, fidelity is 1. As the two cavities move apart, the spectra of the two cavities no 
longer overlap. Thus, there is a small probability of photon detection at d 2 when the dipoles 
are in the state \mm). This results in a loss of fidelity. Surprisingly, however, the fidelity 
does not continue to decrease, but instead increases back to 1 at some value of Aui s . 

As we keep increasing Auj s further, for a certain value of the laser frequency oo, both rf 
and rf are 0. If a detection event is observed in detector d 2 , then the state of the two QDs 
collapses to 

\*f)di P oies = ^[(ar? - /3r?)\mm) + ar™\mg) - f3r™\gm)] (13) 

where N 2 = \ar™ — (3r™\ 2 + |ar'5™| 2 + \/3r™\ 2 . In this special case there is a sec- 
ond matching condition, given by ar™ = flr™, that again projects the two dipoles onto 
|\]/_) = (\gm) — \mg))/y/2. It is this second matching condition that results in the fidelity 
of 1 at the second peak. Our optimization algorithm naturally detects these two optimal 
regions, and gives us the best performance in the intermediate regime. Thus, given any 
set of operating conditions we have the ability to determine the best set of amplitudes and 
input frequencies. We note that in many cases fidelities exceeding 0.95 can be achieved 
even with an 60 GHz detuning, which is more than half a cavity linewidth. The fabrication 
of cavities with resonance frequencies that are repeatable within a linewidth is well within 
current technological capabilities. 




Aw s /k 

FIG. 5: a) Optimized fidelity as a function of Aw s for different values of Aj. A 2 = 

We can also consider what happens when we have both cavity detuning and dipole detun- 
ing. In Fig. 6, we plot optimized fidelity as a function of cavity detuning Au s , and dipole 
detuning Ai. A maximum fidelity of 1 is obtained when Au s = 0, which represents the case 
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FIG. 6: (a) Fidelity as a function of cavity separations Aw s and dipole detuning Ai 

when the two cavities have the same resonant frequency. For Au s < g (where g=20GHz 
in the plot) the fidelity is largely independent of the detuning A x and is only determined 
by cavity separation. When Au s becomes larger, fidelity increases again due to the second 
matching condition and a dependence on the dipole detuning now becomes apparent. This 
dependence on detuning comes about from the fact that the second matching condition is a 
function of Ai, as illustrated in Fig. 5. From Fig 6, fidelity is over 0.75 for a cavity linewidth 
separation(100 GHz) of the cavities even over a wide range of dipole detunings. Thus, we 
can use the protocol to obtain high fidelities even if the cavities and dipoles are detuned. 

VI. 6. EXCITON-BIEXCITON IMPLEMENTATION 

The protocol to generate entanglement between two QDs relies on considering the QD as 
a three level system. One of the ways to implement this three level system is by making use 
of the excitonic and biexcitonic transitions. A QD consists of three states: the ground state, 
an exciton state X consisting of a single electron-hole pair within the QD and a biexciton 
state XX which is formed when two electron-hole pairs are trapped inside the QD. The 
recombination of an electron-hole pair in the XX state generates the biexciton XX photon. 
Similarly, the recombination of an electron-hole pair in X state generates the X photon. The 
X and XX photons have different energies due to the coulomb and exchange interactions 
between the carriers. The typical energy separation between the two lines is lmev[36]. Thus, 
we can make use of this difference in energies to spectrally isolate the two lines. 

The schematic of the QD as a three level system is shown in Fig[7j We identify the three 
states of the QD as the ground, X and XX states. We are free to assign these three states 
as \g), \m) and |e) in a variety of different combinations. In fact, there are several ways 
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to assign these levels, but probably the most convenient approach is given in the inset of 
Fig |7| In the fi gure we have identified the ground state of the QD as state |m), the single 
exciton state as state \g), and the bi-exciton state as state |e). This choice of the level 
configuration has a number of advantages. First, single qubit operations between \g) and 
\m) can be directly applied by pulses resonant with the single exciton transition. Second, by 
placing the bi-exciton transition on resonance with the cavity, we can enhance the exciton 
to bi-exciton transition to get DIT, while at the same time suppressing the single exciton 
lifetime in order to increase the coherence time of the qubit. This is illustrated in Fig [7} 
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FIG. 7: QD as a three level system 



We assume the biexciton transition to be on resonance with the cavity frequency. This 
is indicated in Fig [7] where the XX transition is in the middle of the cavity spectrum. The 
X transition line is detuned from the cavity by 5x- The vacuum Rabi frequencies of X and 
XX transitions are given by gx and gxx respectively. Similarly, the decay rates of the two 
transitions is given by 7x and ^xx respectively. 

In the exciton bi-exciton scheme the degree of cavity enhancement directly impacts our 
ability to create an entangled state. This is because both the exciton and bi-exciton are 
strongly radiative states, and the only way to enhance one while suppressing the other is 
to use cavity lifetime modification. In other qubit implementations, such as dark state 
excitons[26j, this is not as much of a problem because selection rules make the dark exciton 
long lived regardless of cavity. To quantitatively address this issue, we first calculate the 
coherence time of the exciton state which is given by solving for the decay rate of cx_ in Eq. 
1. The coherence time of the qubit is given by 

1 

T 2 



r 



X 



9x K 



X + K 



(14) 



where we have added the dipole dephasing rate 1/T 2 to the decay rate. From the above 
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equation one can see that increasing 5x decreases the decoherence rate until it finally satu- 
rates at a minimum value of 7x + 1/T 2 . At this point, increasing the detuning of the exciton 
will not help as we are limited by non-radiative and dephasing processes. 

The coherence time of the dipole should be compared to an appropriate time scale in 
order to determine if entanglement can be generated. Although there are a number of 
different factors that should be considered in this comparison, the minimum requirement 
for generating entanglement is that the duration of the entangling pulse should be shorter 
than the coherence time of the qubit. If this is not the case, the qubit will begin to decohere 
before the entangling pulse has finished interacting with the cavity-dipole system, and there 
is no hope of generating high-fidelity entanglement. In previous work in Ref [17J, it has been 
shown that when the pulse is resonant with the dipole, it must be much longer than the 
modified spontaneous emission lifetime of the dipole in order to be monochromatic. Thus, 
in the worst case when the dipole is resonant with the cavity we need l/r p << g"xxl K - We 
thus argue that an important figure of merit is the ratio of the coherence time of the qubit 
to the entanglement pulse width, given by 

N ent = &± (15) 

KL x 

This ratio determines the maximum number of entanglements that can be performed before 
the system decoheres. If N ent > 1, there is enough time for the pulses to finish their 
interaction with the QDs before the system has decohered. Otherwise, the QDs will start 
to decohere before the pulses have finished their interaction and high fidelity entanglement 
will be impossible. 

For calculations, we choose experimental values taken from the paper of Hennessy et. al. 
[57] which investigates the coupling of an Indium Arsenide (InAs) quantum dot coupled to a 
photonic crystal cavity patterned in Gallium Arsenide (GaAs) by electron beam lithography. 
This experimental work reports g=20 GHz and cavity linewidth of 25 GHz which corresponds 
to a Q of 13300. However, the cavity linewidth is the bare cavity Q which corresponds to the 
decay into the leaky modes. In order to achieve critical coupling with the cavity, we need 
another in-plane mode with a decay rate equal to the bare cavity decay rate. This mode can 
be implemented in a photonic crystal as a waveguide coupled to the cavity. Thus, the total 
decay rate of the cavity is double that of the bare cavity decay rate. Hence, we use k = 50 
GHz in our calculations. We use gx = gxx = 9- For values of T 2 we use 2 ns, which are 
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appropriate values for InAs QDs[38]. For these values, g 2 xx l ' K ls 8 GHz , Y x is 0.93 GHz and 
N ent is 8.6. The fact that N ent > 1 ensures that we can complete an entanglement operation 
well before the QDs have decohered. 

For a cavity linewidth of 50 GHz, the exciton line lies outside the cavity spectrum(5x = 
250GHz). However, the exciton line still couples to the cavity and we cannot ignore the 
presence of the extra transition coupled to the cavity. So, we cannot substitute for g as 
in Eq 2 in order to obtain the cavity reflection and transmission equations when the QD is 
in state \m). We need to use the vacuum Rabi frequency as the value for g to obtain the 
values of r™, t™, r™ and t™. The changes in the transmission and reflection coefficients will 
modify the final state of the QDs and hence the fidelity of the system. 

In general we cannot assume that the XX transition is not detuned from the cavity 
spectrum. In order to see how robust the biexciton-exciton protocol is dipole detunings, we 
define the detunings of the XX transition lines from their cavities as Sxxi and Sxx2 • In Fig 
[8] we plot the dependence of fidelity on dipole detunings Sxxi and Sxx2 for the above case. 
For both Sxxi = and Sxx2 = 0, fidelity is 1 as expected. When we increase Sxxi and 
Sxx2, the transmission and reflection coefficients are modified due to the coupling of the X 
transition to the cavity. This lowers the fidelity of the output state. The drop is fidelity 
is gradual and for a cavity linewidth separation of the dipoles from the cavity resonance ( 
50 GHz), fidelity drops to only 0.96. As we further increase the detunings to 100 GHz, 
fidelity drops to 0.85. Thus, even for large detunings between the cavities and the dipoles, 
reasonable high fidelity(0.85) states of the QDs can be obtained. Thus, the exciton-biexciton 
scheme can be used to create entanglement between QDs even if the exciton line couples to 
the cavity. The performance of the protocol can be further improved by fabricating cavities 
with high quality factors. 

VII. 7. CONCLUSIONS 

In conclusion, we have shown that one can achieve high fidelity entangled states between 
two dipoles, even when their emission frequencies are different. The method is robust to 
dipole and cavity frequency mismatch. Efficiency loss for a cavity linewidth change in 
dipole detuning is about 50% for a constant fidelity. Therefore, relatively high fidelity can 
be obtained over a large range of dipole detunings without significant loss of efficiency. The 
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FIG. 8: Fidelity as a function of dipole detunings 5i and 62 for the exciton-biexciton model of a 
QD 

development of protocols that are robust to these imperfections is extremely important for 
semiconductor based implementations of quantum networks. 

The authors acknowledge the support of Army Research Office under the Grant No. 
W911NF0710427. 



VIII. APPENDIX 

The input coherent fields are \a) and \f3). Both the dipoles are initialized in a superpo- 
sition of states \g) and |m). Thus, the initial state of the system is 

\%) = \a)\P){\9) + H){\g) + \m))/2 (16) 

The coherent states can be replaced by their corresponding displacement operators to 
account for all order of photon numbers. Thus, 



|^) = i e ^- Q * a -)|0) a!?i e^-^ e -)|0) Ci „ 

2 (17) 

(\gg) + \g m ) + \ m g) + \mm)) 

The input fields after interactions with the cavity-dipole systems are transformed accord- 
ing to Eq 2. Thus, when the dipoles are in state \gg) 



r 9f.\ 
1 2^ out 



2 out 



(18) 
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Similar transformation equations apply when the dipoles are in the states \gm), \mg) 
and | mm). The reflected field from the two cavities is mixed on a 50/50 beamsplitter that 
applies the transformation: 



aL-(d 1+ d 2 )/V2 (ig) 
cL - (di - d 2 )/^2 

Applying the cavity and beamsplitter transformations on the initial state |\I/j) we get 



\% 9 ) -D(^l^ )D{ ^±^ )D{a tl )D (^ 2 ) 

v2 V2 (20) 

out )d ut \gg) 

\^gg) is the state of the output modes for the dipoles in state \gg). This state can be 
split up into the modes of detector di, h out and d out and detector d 2 . Thus, 



l*sa> = \^oo)\^bb)\99) 

l« = D ( ar ? + ^ )jD «)D(^)|0) dl ^ t , dotit (21) 



V2 



\ipgg) is the state of the output modes at detector di and the transmitted modes b out and 
dout when the dipoles are in state \gg). \n gg ) is the field amplitude at detector d 2 when the 
dipoles are in state \gg). Similarly, we can obtain the field amplitudes \^ gm ), \^ mg ) and 
\^mm) when the dipoles are in states \gm), \mg) and |mm). The final state of the system is 
given by 

\Vf) = \V gg ) + |^ Sm ) + \^ mg ) + |* mm ) (22) 

These states \^ gm ), \^ mg ) and \^ mm ) can be further decomposed on similar lines to Eq 
21 to obtain the field amplitudes \ip gm ) and \fi g m)M>mg) and |/i mg ) and \if) mm ) and |/i mTO ) 
respectively. 

We define the projection matrix M as Y^n=i \ n )d 2 ( n \- M can also be written as 
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n=0 

= i-\o) d2 (o\ 



(23) 



tr iBMB) {{M\* f )(9 f \M)} 
PdWOleS tr (dipoles&fields) { (M | | M) } [ ^ } 



,_ ( , r | t n&Ma) {(M\^ f )(^ f \M)} (25) 
'*r (dipoles&fields) { (M| (* / 1 M) } 

The denominator is the probability of getting a detection at detector d 2 . We identify this 
as efficiency 77. 



F= (^-\p dipo ies\^-) 

= ^ | tr( fi eids){|^/)(^/|}-tr (ficlds) {(0|M/ / )(M/ / |0) d2 } |iif 

F\ — F 2 



The individual terms can be evaluated to give 



4 2 2 

F 2 =i e -l^l 2 + ie-l^-l 2 - \e-^\ 2+ \^l\iP mg \^ gm ) 

o o 2 

- ^e-^-l 2+ l^l 2 )/ 2 (^ m |^ s > 

^ — _ jg~M 9 g _|_ g — Cjm _|_ g — Cmj _|_ g~Mmmj 

Thus, the complete expression for fidelity and efficiency can be obtained. 
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